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ABSTRACT 

We present new ground-based data following up on the HST discovery of low-redshift 
Lyman a absorption in the sight-line to the quasar 3C273. Our goal is to investigate 
the relationship between the low-column-density absorbers and higher column- density 
objects such as galaxies or H II regions. Narrow-band filter observations with a corono- 
graph show that there are no H II regions or other strong Ha line-emitting gas within a 
12 kpc radius of the line-of-sight to the quasar, at the velocities of three of the absorbers. 
Broad-band imaging in Gunn r shows that there are no dwarf galaxies at Virgo distances 
with absolute magnitude above Mb ~-13.5 and within a radius of 40 kpc from the line- 
of-sight to the quasar. Finally, we present fiber spectroscopy of a complete sample of 
galaxies within a radius of 1°, down to an apparent magnitude of B^19. Analysis of this 
sample, combined with galaxies within 10 Mpc of the quasar line-of-sight taken from 
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the literature, shows that the absorbers are definitely not distributed at random with 
respect to the galaxies, but also that the absorber-galaxy correlation function is not as 
strong as the galaxy-galaxy correlation function on large scales. We show that our data 
are consistent with the hypothesis that all galaxies more luminous than 1/10 L* have 
effective cross-sections (for association with absorbers whose neutral-hydrogen column- 
density (Log(NH)) is >13.0), of between 0.5 and 1 Mpc. We also show a clear case 
of a Lyman a absorber which has no galaxy brighter than M B =-18 within a projected 
distance of 4.8 Mpc, and discuss the possibility that Lyman a absorbers are destroyed 
in a rich galaxy environment. 

Subject headings: cosmology — interstellar:matter — galaxies:redshifts 
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Introduction 



Understanding the origin and evolution of struc- 



turp in the universe remains one of the most, funda- 
mental and active challenges of current astrophysical 
research. As the evidence in favor of a cosmologi- 
cal origin for the narrow, displaced absorption lines 
in QSO spectra became overwhelming, it also became 
clear that both the metal-line systems and the Lyman 
a systems are invaluable tools for the study of some 
aspects of this problem. Since ground-based Lyman a 
studies refer only to redshifts ;>1.6, they complement 
studies of galaxy clustering properties, the majority 
of which involve redshifts much less than this. How- 
ever, precisely because the redshift regimes have been 
so different and because it has not been at all clear 
what relation exists between the typical low-column- 
density Lyman a absorbers and galaxies, these two 
approaches have remained disjoint. It was somewhat 
unexpected, but pleasing, that low-redshift Lyman a 
absorbers were found in sufficient numbers to enable 
meaningful studies of the evolution of the Lyman a 



absorbers and their relation to galaxies (Morris et al 



1991, Bahcall et al. 1991b). This has presented as- 



tronomers with the opportunity to join these two lines 
of investigation. 

There are two levels at which such such attempts 
can be carried out: 1) Purely statistical investiga- 
tions aimed at comparing the clustering properties 
of galaxies and Lyman a absorbers, and 2) Investi- 
gation of individual cases in which the possibility of 
establishing the presence (or absence) of a clear link 
between a Lyman a absorption line and something 
we could call a "galaxy" presents itself. Preliminary 
discussions along these lines may be found in papers 
by iBahcall et al. 1992a|, |Bahcall et al. 1992b| and 



by palzer 1992 . The present paper is a first attempt 
to pursue both these approaches along the sightlinc 
to 3C273. The remainder of this paper is organized 
as follows: In § ^] we describe the different sets of 
observations we have assembled to investigate the en- 
vironment of the Lyman a absorbers along the 3C273 
sightline. In § [3] we analyze them for possible associ- 
ations or lack of associations of individual Lyman a 
absorbers with galaxies, and also give some statisti- 
cal analysis of the clustering properties of the Lyman 
a absorbers with galaxies. In § [T] we discuss these 
results in light of current models for the Lyman a ab- 
sorbers and provide a brief summary and suggestions 
for further work. 



Throughout this paper H is taken to be 80 km/s/Mpc, 
the distance to the Virgo cluster is taken to be 16.0 
Mpc (a distance modulus of (m-M)vi rg o=31.02) ( |Ja-| 
coby et al. 1992), and qo is taken to be 0. 



2. Observations and Reduction 

It has long been realized that imaging of the gas 
directly responsible for the low-column-density Ly- 
man a absorbers is well beyond the reach of current 
technology. Specifically, the neutral-hydrogen column 
densities of order 10 13 - 10 14 cm" 2 that GHRS de- 
tected toward 3C273 are about 4 or 5 orders of mag- 
nitude below what can be imaged in 21 cm emission, 
even without taking into account the powerful radio 
background contributed by 3C273 itself. The Ha re- 
combination surface brightness associated with these 
neutral-hydrogen column-densities is also several or- 
ders of magnitude below what is feasible to detect, 
unless the incident flux of ionizing photons is several 
orders of magnitude higher than that expected from 
the integrated background radiation. 

However, it has frequently be suggested that the 
Lyman a absorbers are intimately connected with, or 
are actual extensions of, entities which can be imaged 
either by means of Ha emission, starlight, or 21 cm 
emission - e.g. dwarf galaxies (Fransson and Epstein 
1982|) or shells of expanding gas (Chcrnomordic anc 



Ozcrnoy 1983 ) or the outer regions of galactic disks 
( Maloncy 1992 ). In the case of dwarf irregulars, for 
example, a very small episode of recent star formation 
might betray the presence of a dwarf irregular whose 
outer envelope produces the Lyman a absorbers. Al- 
ternatively, expanding shells of gas might produce Ha 
emission via collisional ionization at a shock front. 

In addition, of course, possible association of indi- 
vidual Lyman a absorbers with specific galaxies, as 
well as statistical studies of absorber-galaxy correla- 
tion can be carried out with a sample of redshifts for 
galaxies in the field surrounding 3C273. 

In this section we describe three such new sets of 
observa tion s of a region centered on 3C273. These 
are: § (2.lJ Coronagraph observations with narrow- 
band filters of a 5' diameter region, § 2.2. Deep broad- 
band imaging of a 17' diameter region, and § 2.3 
Fiber spectroscopy of a 2.2x1.6° region down to a 
limiting magnitude of B=19. We describe the analysis 
of these three sets of observations in § ||] 
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2.1. Coronagraph Observations with Narrow 
Band Filters 

Observations of a 5f3x5.'3 region (radius «12 kpc 
at Virgo) around 3C273 were obtained during 1992 
Februa ry 3-7, with the Unive rsity of Hawaii Corona- 
graph ( Vilas and Smith 1987 ) on the Las Campanas 
2.5m duPont Telescope. A thinned 1024x1024 Tek- 
tronics CCD was used, binned 2x2, giving a scale of 
1.23"/pixel. The coronagraph blocking mask had a 
diameter of 5". Data were obtained with a Gunn r fil- 
ter, and also 5 specially acquired filters, 3 with width 
13.5A, centered at 6586. 2A, 6598. OA and 6756.0A 
(hereafter referred to as VN1, VN2 and HN), and 2 
with width 25A, centered at 6643.2A and 6718.9A 
(hereafter referred to as VB and HB). (The above 
widths and central wavelengths are quoted for an f7.5 
beam and a temperature of 15° C). The narrow-band 
wavelengths were chosen to match the redshifted po- 
sition of Ha at the velocities of the 2 Lyman a ab- 



sorbers listed in Morris et al. 1991 at velocities corrc- 



sp onding to the Virgo cluster ( Bingclli, Sandage and 
Ta mmann 1985 ), and one absorber at 1251A. This is 
the lowest redshift, strong Lyman a system beyond 
the Virgo cluster. The observing procedure involved 
cycling through the 6 different filters, with exposure 
times of 10 mins for the r and 25A filters, and 20 
mins for the 13. 5 A filters. Seven such cycles were 
completed over the 5 night run. During the observing 
run, it was discovered that the narrow-band filter HN 
had slightly non-parallel faces, resulting in detectable 
'ghost' images offset from bright stars and also 3C273. 
In an attempt to minimise the effect of these, this fil- 
ter was rotated though 90° between each night. 

For calibration purposes, observations were also 
obtained of Mkn 49 (an emission line galaxy in the 
Virgo cluster with radial velocity 1524 km/s, and 
hence with Ha line-emission within 2A of the peak 
of VN2) and M87, and also a number of bright stan- 
dard stars. 

The images were reduced using IRAF^. The reduc- 
tion steps were: bias subtraction, division by a flat 
field taken on the same night as the observations, ro- 
tation and shifting of the images to match a reference 
image, sky subtraction and averaging together of im- 
ages taken with the same filter. It was found after the 
run that refocusing the coronagraph between taking 
the flat fields and the data meant that the flat field di- 



3 IRAF is distributed by NOAO, which is operated by AURA 
Inc. under contract to the NSF 



vision left significant features in the data, both at the 
edges of the coronagraph field and also throughout 
the data at the location of what are presumed to be 
dust particles on the coronagraph optics or CCD win- 
dow. This problem was particularly noticeable for the 
data taken on the last night of the run. However, the 
residuals are greatly reduced in the combined data. 

Continuum sources were removed from the narrow- 
band images by subtracting off a scaled version of 
the two 25A filter observations. We investigated scal- 
ing methods, including measuring stars or galaxies in 
the images, but found that the scaling derived was 
consistent with simply subtracting off the average of 
the 25A filters from each 13. 5A observation. That 
is, there was no evidence for a significant continuum 
slope or calibration difference across the 150A region 
of interest, and the factor 2 shorter broad-band expo- 
sure fairly accurately balanced the higher throughput 
of the 25 A filters. This apparent consistency may 
be fortuitous: due to the small field (small number of 
galaxies), and the errors in measuring the flux of faint 
galaxies. Figure [i] shows the sum of the 25 A data, and 
the continuum subtracted 13. 5A data for 3C273. As 
can be seen, the PSF is a function of angle off-axis, 
and becomes quite broad and asymmetric at the edge 
of the field, due to aberrations in the coronagraph 
optics. 

2.2. Wide Field Imaging with COSMIC 

A mosaic of Gunn r band images of the field around 
3C273 was obtained on 1992 February 23, at the 
Palomar 5m, with the prime focus COSMIC sys- 
tem. This recently commissioned camera contains 
a thinned 2048x2048 Tektronics CCD. Due to poor 
seeing, this was read out with 2x2 binning giving a 
scale of 0.56"/pixel. Exposures were taken roughly 
centered on 4 positions offset from 3C273 by 5.'4 NE, 
NW, SE and SW. Four exposures of 5 minutes each 
were obtained at each position. During the observa- 
tions, the moon rose, and so the background level of 
the images varies by almost a factor 2. 

The data were reduced using IRAF. After bias 
subtraction, the data were flat fielded using images 
of the dome. This left significant large scale struc- 
ture in the images, and so a 'skyflat' was constructed 
from the combined data images. This flat was me- 
dian smoothed to remove any small scale structure. 
Because of a region of very low sensitivity near the 
center of the CCD, and also the location of a very 
bright star coincidentally at the same place in one set 
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of images, the flat fielded data still shows a weak neg- 
ative feature near the center of each image. The sky 
background was determined for each image separately 
and subtracted. The offsets between the images were 
determined by measuring the positions of the QSO 
and 2 bright nearby stars, which were common to all 
images and the data were then mosaiced together, 
giving a combined image with diameter 17.'2 centered 
on 3C273. The resulting image is show as figure 0. 

2.3. Fiber Spectroscopy 

During f992 February 8-fO, spectra of objects in a 
2.2 xf. 6° rectangle surrounding 3C273 were obtained 
with the Fiber Spectrograph at the Las Campanas 
2.5m duPont Telescope. ( [Bhectman 1992| ). This sys- 
tem has 128 fibers which are manually plugged into 
an aluminum plate over a 1.5x1.5° field. The fibers 
have a projected diameter of 3'.'5. They feed the slit 
of a floor- mounted spectrograph. A 600 1/mm grating 
was used, with the 2D-Frutti photon-counting detec- 
tor, giving a resolution of 8.6A FWHM. Three fields 
were observed for about 6000 seconds each, offset EW 
from each other. 

Objects observed were chosen from a database pro- 
duced by scanning an UKST IllaJ plate of the region 
with the APM scanning machine at Cambridge. The 
APM produces a catalog of all the objects on a plate, 
with estimates of isophotal magnitudes, size and a 
'sigma' parameter that measures how much the image 
parameters differ from those of stars with compara- 
ble magnitude. We chose objects to observe from this 
catalog with 'sigma'>3.0. This includes many fairly 
compact objects, and resulted in a rather high con- 
tamination by stars, but also means that we found 
a number of compact galaxies that would otherwise 
have been missed. For each 1.5x1.5° fiber field, a 
magnitude sorted list of candidate galaxies was pro- 
duced. Due to restrictions on the minimum fiber sep- 
aration, 330 out of a possible 336 object fibers were 
used, with 16 fibers per field set on blank sky. 

The fiber spectra were extracted and reduced using 
the IRAF apextract package. The spectra were first 
traced and extracted. Then fiber-to-fiber throughput 
differences were corrected with a flat field image, and 
the fibers were wavelength calibrated and rebinned to 
a common linear wavelength scale. The wavelength 
calibration used an arc spectrum to determine the 
non-linear relationship between wavelength and pixel 
number, after which a zero point shift was measured 
for each fiber using the strong sky lines in the data. 



Finally, the sky emission was subtracted from each 
fiber using an unsealed template constructed from the 
16 sky fibers in each frame. 

Classification and radial velocity measurements for 
each spectrum were done in two ways. First each 
spectrum was inspected by eye, and classified as ei- 
ther a star, galaxy or unknown. Then a radial veloc- 
ity was determined by measuring either the position 
of the 4000A break, or that of the [Oil] A3727.6 fea- 
ture (actually a doublet, but unresolved in our data). 
A subjective assessment was also made of the reli- 
ability of the resulting radial velocity, dividing the 
sample into 'possible' and 'definite'. All the spectra 
were also analysed using the fxcor routine in IRAF. 
Each object was cross-correlated with three different 
templates: (a) a template made up from the 27 best 
S/N stars in the data, (b) a template made up from 3 
high S/N late-type stars, and (c) a template made up 
from 16 emission line galaxies (all shifted to their rest 
frames) in the data. For each spectrum, the correla- 
tion with the highest peak was then selected. A re- 
assuringly close match was found between the by-eye 
classification and velocities and the results from fxcor. 
The resulting histogram of velocities was inspected. 
Apart from 4 low S/N or hot stars (for which a good 
template was not constructed), the stellar velocities 
found by fxcor had a distribution well represented by 
a gaussian with zero mean and an approximate dis- 
persion (a) of 85 km/s. We take this to be a reason- 
able estimate for our radial velocity uncertainties. It 
was also found that the subjective 'possible' category 
in the by-eye radial velocity measurements matched 
rather well to a cross-correlation peak height less than 
0.3 returned by fxcor. Apart from one outlier, the 
difference between the fxcor velocity and that mea- 
sured by eye for the galaxy identifications with cross- 
correlation peak heights above 0.3 also was fairly well 
fit by a gaussian with a of 85 km/s. 

In the end we obtained 129 definite galaxies, 43 
possible galaxies, 86 definite stars, 4 possible stars, 
10 fibers that had to be unplugged due to over- 
illumination (which were hence almost certainly stars), 
37 fibers which showed no flux (either very low surface 
brightness galaxies or positional errors) and 21 spec- 
tra which showed flux but which were unclassifiable 
either as stars or galaxies (based on very low cross- 
correlation peak heights from fxcor and visual inspec- 
tion). We should reiterate that a much higher 'suc- 
cess' rate in finding galaxies could have been achieved 
by raising the cutoff 'sigma' value used in choosing 
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galaxies from the APM scans, at the cost of missing 
some compact galaxies. 

We present in table [l] the resulting galaxy red- 
shifts. For each object, the table lists the RA and 
Dec, an approximate B magnitude and the heliocen- 
tric radial velocity. The magnitudes were calculated 
using the APM isophotal magnitudes measured from 
the plates, crudely calibrated using the B magnitudes 
listed by Stockton 1980 and Salzcr 199^ for objects 
in the scanned region. They could be in error by as 
much as 0.5. 

A number of the brightest galaxies in the field were 
not included in the fiber survey. Those with known 
redshifts within the survey region were added to the 
sample (4 objects taken from the May 5 1990 version 
of the CfA Redshift Catalog, |Huchra 1990|). This 



gives a total sample of 176 galaxies with redshifts 
within the survey region which is roughly complete 
to a B magnitude of 19.0. 

We plot the results of the survey in a number of 
projections. Figure p| shows a redshift histogram of 
the galaxies, figure [|shows the distribution of galax- 
ies on the sky, while figure || shows the pie-diagrams 
obt ained , 



3. Analysis 

In this section we go through each of the three 
observations described in § ^\ in t urn, deriving con- 



straints on the absorbers. In § 3.1. we derive limits on 
Ha line emission at the velocities of the narrow-band 
filters and in § 3.2. we calculate the maximum abso- 



lute magnitude a dwarf galaxy could have and remain 
undetected in our broad-band imaging. A long de- 
scription of the correlation analysis between the Lya 
absorbers and the galaxies found in the fiber survey is 
given in § 3.3. , in which the various available absorber 



and galaxy sub-samples are discussed, and two al- 
ternative extreme hypotheses for the absorber -gala xy 



correlation function are tested. Finally in § 3.4.| we 
discuss some particular aspects of the absorber-galaxy 
distribution found in the 3C273 sightlinc. 

3.1. Flux Limits for Ha Line Emission 

The final continuum-subtracted coronagraph im- 
ages have a measured rms dispersion (away from 
residuals due to bright stars) of 0.0018 DN/pixel/s. 
From calibration observations of HD84937, this is 
equivalent to a la flux limit of 2. 10 — 18 ergs/cm 2 /s/D", 
equivalent to an emission measure of approximately 



2.8 pc cm~ 6 . By blinking the images, it can be seen 
that none of the objects visible in the broad-band 
images have emission line fluxes > 3c. For compari- 
son, the VN2 image of Mkn49 showed a peak Ha line 
flux of 3.10~ 14 ergs/cm 2 /s/D", i.e. 5,000 times higher 
than our 3a detection limit. One can also perform the 
following thought experiment. What would the Orion 
nebula look like if placed at the distance of the Virgo 
cluster? As discussed in Kennicutt 1984, the Orion 



nebula is actually a relatively low luminosity H II re- 
gion, with an Ha luminosity of only 10 37 ergs/s. Nev- 
ertheless, if placed at a distance of 16 Mpc, it would 
still have a flux of 1.3 10~ 15 ergs/cm 2 /s (and would 
be unresolved - the nebula has a diameter of 5 pc, 
while at Virgo, 1" corresponds to approximately 80 
pc). Thus it would be a factor 220 brighter than our 
3a limit. One can also calculate that the Stromgren 
sphere around a single main sequence star of spec- 
tral class ^Bl/D" would be detectable at the 3a level 



( |Allen 1973| ). 

Unfortunately, the expected surface brightness of 
an optically thick slab of hydrogen, simply bathed 
in the local UV background would not be d etectable. 
Taking the limit on the UV background from Bongaila"] 



Bryant and Cowie 1989, and using the formulae from 



Ostcrbrock 198S , one finds a surface brightness of 
< 3.10" 19 ergs/cm 2 /s/D", about a factor 20 below 
our 3a limit. 

Higher spectral and spatial resolution data have 
been taken by T. Williams (private communication) 
using the Rutgers Fabry-Perot system at the CTIO 
4m telescope. This data has not been fully analysed, 
but should produce even lower surface brightness lim- 
its (or detection) than our data. 

3.2. Limiting Magnitude for detection of Low 
Surface Brightness Dwarf Galaxies 

The main motivation for taking the COSMIC im- 
ages was to determine whether there were any dwarf 
galaxies within a 40 kpc radius from the line of sight 
to 3C273 which are too faint to be visible on the POSS 
sky survey plates. Examples of low surface brightness 



dwarfs in the Virgo cluster can be seen in Sandage 
and Bingclli 1984. In order to determine how faint 



a dwarf galaxy, which had morphological properties 
typical of those found in the Virgo cluster, would be 
detectable in our image, we used the IRAF artdata 
package to insert artificial dwarf galaxies into the data 
array. The magnitude scale was derived by assuming 
that the sky background in the Gunn r band (before 
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the moon rose) was 21.5 (Massey 1990). B absolute 
magnitudes were converted to Gunn r assuming a B-R 
color of 1.0 (B-V of 0.7), and a conversion from R to 
Gunn r of r=R+0.43+0.15(B-V) ( [Kent 1985| ). Thus a 
dwarf galaxy with Mb=-15.5 was taken to have M r =- 
16.0. The Virgo distance modulus was taken to be 
31.02. Dwarf galaxy properties were taken from Bin 



gel i, Sandage and Tammann 1985. In particular, typ- 



ical exponential scale lengths^ for Virgo dwarfs were 
taken to be 2-4 kpc. Figure oshows the same data as 
figure but with three dwarf galaxies added. One to 
the NE with Mb=-14.5 and scale length 4.4 kpc, one 
to the NW with M B =-13.5 and scale length 2.2 kpc, 
and one to the SW with Mb=-13.5 and scale length 
4.4 kpc. These experiments demonstrated that we 
would be able to detect dwarf galaxies as faint as -13.5 
at the distance of Virgo, having morphological prop- 
erties similar to those found near the cluster center. 
For comparison, the dwarf galaxy illustrated in figure 
2, panel 4, labelled T5°47' of |Sandagc and Bingelli 
198^ has an absolute B magnitude of -14.8. 

A complementary approach is to search for H I 
21 cm emission from Virgo dwarf galaxies. Re- 
cently, van Gorkom and her collaborators have set 
extremely low limits with the VLA over a 40x40 
arcmin field centered on 3C273 and over a velocity 
range of about 1000 km/s centered on 1300 km/s to 
a 1 sigma column-density limit of approximately 10 19 
cm~ 2 ( van Gorkom 199*3] , |van Gorkom et al. 1993 ). 



3.3. Correlations between Lyman a Absorbers 
and Galaxies 

We have shown there are no H II regions, or other 
strong Ha line-emitting gas, or dwarf galaxies near 
the Virgo absorbers. Having determined the rcdshift 
distribution of galaxies near the sightline to 3C273, 



(§ 2.3.) we would now like to address the statisti- 
cal question of the degree to which the Lyman a ab- 
sorbers are correlated with luminous galaxies. If they 
are correlated, is the Lyman a absorber-galaxy cor- 
relation the same as the galaxy-galaxy correlation? 
A cursory inspection of figure ^| is enough to show 
that there will not be a simple answer to this ques- 
tion. While there do seem to be Lyman a absorbers 
associated with clumps of galaxies (e.g. the Virgo 
absorbers, or the set of absorbers around z=0.02- 
0.03), there are also absorbers in conspicuous 'voids' 
(at z=0.06-0.07), and there are no absorbers associ- 



ated with the prominent excess of galaxies at z=0.078. 
We consider statistical tests for various assumptions 
about the correlation between the Lyman a absorbers 
and the galaxies near the sightline towards 3C273 for 
which we have redshift information. Ultimately, the 
goal should be a quantitative and complete statistical 
description of the clustering properties of the Lyman 
a absorbers themselves and their correlation with var- 
ious types of galaxies, clusters, voids etc. Given both 
our rather meagre understanding of this problem and 
the small data set, we shall concentrate on testing 
the following two extreme null hypotheses about the 
Lyman a absorber clustering properties: 

1. The Lyman a absorbers are uncorrelated with 
galaxies and are randomly distributed. (The 
second part of this assumption necessarily im- 
plies the first, but the converse is not necessarily 
true: The absorbers could be correlated among 
themselves but be uncorrelated with galaxies) 

2. The Lyman a absorbers are correlated with 
galaxies in the same way that galaxies are corre- 
lated. More precise formulations of this hypoth- 
esis depend upon the particular test applied, as 
described below. 

In carrying out most of the tests described below, 
it is necessary to compute the 3-dimensional distance 
between every absorber-galaxy or every galaxy-galaxy 
pair. The question arises as to how to compute the 
component along the line of sight, since departures 
from a perfectly smooth Hubble flow distort the map- 
ping of redshift onto radial distance. With no infor- 
mation other than the angular coordinates and red- 
shifts for the objects, we cannot uniquely determine 
the separation along the line of sight for any individ- 
ual pair of objects. For purposes of statistical tests 
we therefore make two different assumptions about 
this component; the degree to which we do or do not 
derive similar results will provide some indication of 
the sensitivity of the test to the uncertainty in the 
estimation of this component: 

1. We simply ignore any departures from Hubble 
flow 



2. We a dopt the formalism of Davis and Peebles 
1983| to estimate this component. These au- 



*Scale Length=Ro, with Intensity oc exp(— 1.6783 X R/Rq) 



thors show how, knowing the two-point correla- 
tion function for the projected distance between 
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pairs, in principal one can invert the integral 
equation relating the projected distance corre- 
lation function to the three-dimensional spatial 
correlation function. However, for our limited 
and rather noisy sample, this is not a very sat- 
isfactory procedure. Moreover, we would like 
statistical estimates for the 3-dimensional sepa- 
ration for each pair for the purposes of carrying 
out other types of tests (e.g. nearest neighbor 
tests). We therefore use the Davis-Peebles for- 
malism as follows: Adopting the assumptions 
described in their paper, the integrand in their 
equation (22) represents the probability that a 
pair with projected separation r p Mpc and ve- 
locity difference tt km/s has a separation in the 
radial direction of y Mpc. We further assume 
their functional form for h(r) in their equation 
(23) , with the parameter F = 1, and adopt their 
functional form for f(V), the probability of the 
relative velocity difference, as well as their ex- 
pression for (j, the dispersion in f(V), i.e. their 
equation (32). For r p » r a or n » H a X r a , 
where r Q is the characteristic correlation length, 
the probability is strongly peaked at y ~ H x n, 
i.e., the pair separation is very likely to be that 
given by assuming a pure Hubble flow. How- 
ever, as both of these inequalities fail to be sat- 
isfied, a second maximum in the probability dis- 
tribution arises at y w whose height depends 
upon the strength of the correlation - i.e., a pair 
with a moderate velocity separation and small 
projected separation may be separated by their 
Hubble flow distance, but if the correlation be- 
tween such pairs is strong, this moderate ve- 
locity separation is more likely to arise from a 
pair at about the same distance from us, with 
strong gravitational interaction between them. 
As noted at the outset of this discussion there 
is no way to unambiguously determine the rel- 
ative separation of any pair along the line of 
sight, but since we are interested in statistical 
applications we adopt the expectation value of 
this probability distribution as our second alter- 
native algorithm. 

A problem with this second approach is that evalu- 
ating the probability distribution for the radial sep- 
aration of the pair requires that we know in advance 
the two-point correlation function between the pair. 
Ideally, we could have dealt with this problem by 
using an iterative approach: for both the absorber- 



galaxy and galaxy-galaxy pairs separately, we start 
with some "fiducial estimates" for the two-point cor- 
relation (e.g. r a — 5.4/i _1 Mpc and 7 = 1.77 ; c.f. 
equation 19 in Davis and Peebles) to compute the ex- 
pectation value of the radial separation for each pair, 
and thus compute the two-point correlation functions. 
With these new, separate best-fit values of 7 and r D 
for the correlation functions for the galaxy-absorber 
pairs and the galaxy-galaxy pairs, we could then re- 
peat the process until the parameters for the two- 
point correlation functions have converged. In fact, 
since the results of our statistical tests do not appear 
to be very sensitive to departures from a pure Hubble 
flow, we have not carried out this iteration, but have 
simply used the single set of parameters (r a — bAh^ 1 
Mpc and 7 = 1.77) in calculating the expectation val- 
ues for both sets of pairs. 

In the following, we refer to these two algorithms 
for estimating the radial separations as the "pure 
Hubble flow" and "perturbed Hubble flow" cases. 
When listing object separations in the following sec- 
tions we will give the separations found from the ex- 
pectation value of the perturbed Hubble flow model 
in brackets following the value for the pure Hubble 
flow model. 

Before carrying out any statistical tests we define 
the two samples and discuss appropriate corrections 
for completeness. 

3.3.1. The Lyman a Sample 

A carefully defined list of Lyman a absorption lines 
is essential to a proper statistical discussion of the cor- 
relation properties of Lyman a absorbers with galax- 
ies. The preferred list would obviously be drawn from 
a homogenous set of observations with the smallest 
detectable equivalent width covering all or most of 
the relevant redshift range. 



Several line lists for the 3C273 sight line have 



been published (Morris et al. 1991 




Bahcall et al. 


1991a 




Bahcall et al. 1991b, 


Brandt et al. 1993 




Bah- 


call et al. 1993 


). These line lists are compared in 



table |2|. However, not only are these lists based upon 
3 different HST spectrograph configurations (GHRS 
G160M, GHRS G140L, FOS G130) but they have 
also been produced by different reduction procedures 
and line-finding and measuring algorithms and with 
differing acceptance criteria for what constitutes a 
"real" line. To investigate the importance of this lat- 
ter source of inhomogeneity we have run the same 
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continuum fitting, absorption line-finding and line- 
measuring software, "JASON", used for the FOS line 
list, which is described in detail by [Bchncidcr et al 



19£|3f| on the data sets of Morris et al. 1991 and 
Brandt et al. 19931. The JASON software was de- 



signed to run on FOS data with an approximately 
gaussian point spread function (PSF). Unfortunately 
this is not a good representation for the GHRS large 
aperture PSF, and also in general the lines were re- 
solved - meaning that the observed line profiles had 
neither the instrumental PSF nor a gaussian shape. 



Thf current version ot JASU1N docs not pcrtorm sucli 



convolutions, and so we nave run the search routines 
assuming a fixed PSF with the correct (non-gaussian) 
shape for the GHRS, but with no account taken for 
resolved lines. This means that the EWs output by 
the JASON software are not accurate, but the detec- 
tion significance levels are approximately correct (see 
Schneider et al. 1993). The data set used is that 



described by Brandt et al. 1993. We tabulate the sig- 
nificance levels from JASON for the Lyman a lines in 
column 5 of table p[ after the positions and EWs pub- 
lished in Brandt et al. 1993. It can be seen that all 



of the 'reliable' lines listed by Morris et al. 1991 arc 
confirmed by the JASON software, with the notable 
exception of the line at A1276.54, which was also not 
found by Brandt et al. 1993. We checked this line 



by running the JASON software on the original data 
used by Morris et al. 1991, obtaining a significance 



lev*l ot 4.5 tor the line. The line finding and fitting 



software used m Morris et al. 1991 was developed by 
R. Carswell and J. Webb, and is described in that 
paper. It used the GHRS PSF convolved with a Voit 
profile with variable width. 

Two other entries in table |^ require special com- 
ment: The li ne at A=1317.08 is ide ntified as Lyman a 
in the list of Bahcall et al. 1991b| whereas in the list 
of Morris et al. 1991 it is identified as Ni II. This is- 



sue has been discussed in detail in Brandt et al. 1993 



wh ) give reasons for preferring the Ni II identifica- 



tion which we adopt. The second case involves the 



line at A=1 393.86. As discussed in detail by pavage 
et i iL 1993 this line appears to be a blend of one of 
the members of the galactic Si IV doublet and an- 
other strong line, whose only plausible indentification 
is Lyman a. The procedure used by Morris et al 



19S(T] in estimating the strength of this line (which in 



5 We thank D. Schneider for kindly making available the most 
recent version of JASON and for instruction in its use. 



volves a detailed comparison of the line profiles of the 
Si IV doublet) is not incorporated into the JASON 
formalism. For this reason we cannot assign a formal 
uncertainty in the line strength. However the resid- 
uals from an unblended fit to the Si IV doublet are 
highly significant. 

As a result of the above considerations, we have 
decided to adopt the following samples of absorbers 
for our statistical tests: (1) For the Lyman a ab- 
sorber "full sample" we adopt the list of 16 Lyman 
a absorbers (and their redshifts) given in Morris e\ 



al. 1991 along with the additional low-redshift line 



(A1224.52) given by Bahcall et al. 1991a; as noted in 
table g, this last line is visible in the GHRS G140L 
spectra but was b elow the significance threshold of 
Morris et al. 1991. This "full sample" is inhomoge- 



nous and/or biased in three senses: i) The high resolu- 
tion GHRS data covers only redshifts above z~0.016; 



below this, only the FOS data of Bahcall et al. 1991b 



and the GHRS G140L low resolution data are avail- 
able. (Observations with the GHRS G160M grating 
in the redshift regime from 0.0 to 0.016 are sched- 
uled for HST Cycle 3). Thus, the full sample may 
be biased against weak low-redshift lines in this red- 
shift range that may be detected by these Cycle 3 
observations, ii) The line of sight toward 3C273 may 
be somewhat atypical in that it passes through the 
southern extension of the Virgo cluster, iii) Some of 
the weakest lines listed as "possible" in the Morris 



et al. 1991 list may not be real. Accordingly, along 
with the full sample of 17 lines, we shall also con- 
sider two subsamples: (2) A "homogeneous sample" 
made of the set of 14 lines found only with the GHRS 
G160M observations using the original Carswell and 
Webb software (i.e. all but the first 3 lines in column 
1 of table ^). (As it happens, this is also equivalent 
to deleting the 3 low-redshift (z^0.016) lines possi- 
bly associated with the Virgo cluster). (3) A "strong 
sample" composed of the set of 10 lines from Morris 



et al. 1991 with -log(P)>7.5 (but including the line 



at z=0. 14658 for which a formal probability estimate 
was not possible due to blending - i.e. the above sam- 
ple with the lines marked 'd' in the comments column 
of table H removed) . We have listed the sample mem- 
bership in the final column of table [|. Note that in 
contrast to the galaxy sample discussed below there is 
no intrinsic observational selection against the higher 
redshift absorbers. 

One could, of course, define further samples. In 
particular, at the request of the referee, we have also 
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run our statistical tests of the complete set of 5 lines 
listed in Bahcall et al. 1993 (i.e. the line list in col- 



uma 8 of table 2). Unfortunately, the number of 



lines in this list is so small that neither of the two 
null hypotheses considered below in connection with 
cloud-galaxy association can be rejected with any sig- 
nificance. For this reason, and in order to keep the 
various combinations of absorption line and galaxy 
samples to managable proportions, we limit our ta- 
bles of statistical results to consideration of the three 
samples defined above. 

3.3.2. The Galaxy Sample 

An appropriate sample of galaxies with which to 
carry out the correlation analysis would be one which 
is complete to some limiting absolute magnitude through- 
out a cylindrical volume centered on the 3C273 sight- 
line (i.e. out to a constant impact parameter) with a 
radius large enough to sample most of the expected 
power in the correlation function and length over (and 
beyond) the full redshift range covered by the Ly- 
man a line sample. The observed sample described 
in § |2.3.| fails this requirement in two obvious respects: 
i) It detects only the more luminous galaxies at the 
higher redshifts, ii) It contains no galaxies with large 
impact parameters at low-redshifts. For some tests 
these deficiencies are probably not important but for 
others they are. Accordingly, we will consider two 
galaxy samples. The first is the sample described in 
which we will refer to as the 'cone' sample. The 



2.3 



second is the cone sample together with all galaxies 
from the May 5 1990 version of the CfA Redshift Cat- 



alog, QHuchra 1990Q within 10 Mpc of the 3C273 line- 
of-sight. This gives a heterogeneous sample with an 
unknown selection function, but is closer to the ideal 
'filled cylinder' than the cone sample. It contains 1498 
galaxies, the vast majority being at Virgo distances, 
and will be referred to as the 'cylinder' sample. 

Thus for each hypothesis tested below there are 6 
galaxy-absorber sample combinations, and 2 possible 
estimates for the distance between every pair of ob- 
jects giving a total of 12 data sets for each statistical 
test. 

3.3.3. Tests of the 1st Null Hypothesis: The Lyman 
a absorbers are uncorrelated with galaxies 

It is obvious that this hypothesis cannot literally be 
true - every sightline that passes close to any galaxy, 
except those utterly devoid of gas, will surely pro- 



duce a detectable Lyman a line, and indeed examples 
of this are already known (Bahcall et al. 1992a, Bah 
call et al. 1992b| ). Nevertheless, in light of the fact 



that the high-redshift Lyman a absorbers show al- 
most no power in their two-point correlation function 



(c.f. Rauch et al. 1992 and references therein) it is of 
interest to see if the present data set does or does not 
exclude this hypothesis, and if it does, how strongly. 

Having formulated this null hypothesis we consider 
two statistics as measures of correlation (or lack of it): 

1. The average over all the Lyman a absorbers of 
the distance between a given Lyman a absorber 
and the N nearest galaxies in the sample, with 
N=l,3 and 5. 

2. The total number of absorber-galaxy pairs within 
a fixed radius R, with R=500 kpc and 10 Mpc. 



In order to see whether the values of these observed 
statistics are such that the null hypothesis can be re- 
jected, we must find the distribution of these same 
statistics for many realizations of the null hypothesis. 
To do this, we carry out 1000 Monte Carlo simula- 
tions in which the same number of absorbers as that 
of the particular Lyman a sample under considera- 
tion are laid down randomly (but follow the "global" 
distribution dN J dz ~ (1+z) 3 as determined by Bah- 
call et al. 1993 ). Redshift limits for the random ab- 



sorbers were 0.03<z<0.151 for the total sample, and 
0.016<z<0.151 for the homogeneous and strong sam- 
ples. 

The results for the nearest-neighbour tests are 
summarised in table ||. We show results for the single 
nearest neighbour galaxy, the mean of the 3 near- 
est and the mean of the 5 nearest. For each case, 
the three columns list the observed mean nearest- 
neighbour(s) distance, the average of the mean nearest- 
neighbour^) distances produced by the Monte Carlo 
simulations, and the number of the 1000 Monte Carlo 
simulations that had a mean nearest-neighbour dis- 
tance less than the observed one. Thus this last col- 
umn, divided by 1000, can be taken as the probability 
that the observed value could arise from a sample of 
absorbers distributed at random with respect to the 
galaxies. 

The most striking result is that there is a less than 
0.1% probability that the average nearest-neighbour 
distance to the single closest galaxy could arise from 
a randomly distributed set of absorbers. This is true 
for all sample combinations, and for either pure or 
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perturbed Hubble flow. For all samples, as one in- 
cludes more galaxies in the nearest-neighbour aver- 
age, the significance drops. After seeing this result, 
we wanted to test whether all the significance came 
from the nearest galaxy, and so ran tests on the sec- 
ond nearest galaxy only (also given in table |^) . One 
can see that the observed mean distance is still signif- 
icantly lower than that expected for a randomly dis- 
tributed set of absorbers, but this could be explained 
by a combination of the highly significant correlation 
with the nearest galaxy and the strong galaxy-galaxy 
two-poin t correlation. This point is discuss ed in some 
detail by Phillips, Disney and Davies 1995 in the con- 
text of bright galaxies found near quasar Mgll ab- 
sorbers. 

We also give in table [| the RA, Dec and redshift 
of the nearest galaxy in the 'cylinder' sample to each 
Lya absorbers of the 'total' sample (assuming pure 
Hubble flow). The final column in this table is the 
minimum absolute magnitude that could have been 
detected in the fiber survey. This shows that for one 
absorber there is no known galaxy with absolute mag- 
nitude above -17.8 within nearly 10 Mpc, (the same 
distance for the perturbed Hubble flow model) and 
that even the nearest absorber-galaxy pair in our sam- 
ple are separated by 350 kpc (240 kpc) . We will return 
to this in 



3.4 



The results for the number of galaxies within a 
fixed radius of each absorber are given in table §. 
These tests are essentially a comparison of the in- 
tegrated two-point correlation function out to the 
given radius (see Mo et al. 1992 for a discussion 
of this point). The numbers given in the table are 
the observed number of absorber-galaxy pairs within 
the given radius, the average of the number of pairs 
found in 1000 Monte Carlo simulations with randomly 
distributed absorbers, and the number of the Monte 
Carlo simulations with a larger number of pairs than 
that observed. Thus the final column divided by 1000 
is the probability that the observed numbers of pairs 
or more would arise from a random distribution of 
absorbers. 

Table |5| shows that there is no significant excess 
of galaxies within volumes of radius 10 Mpc centered 
on the absorbers compared to a random distribution, 
apart from the cylinder/total subsample. For this 
combination there are a large number of pairs between 
the Virgo absorbers and the many Virgo galaxies in 
the CfA catalog. This result may be interpreted as 
saying that it is surprising to find 3 out of 17 absorbers 



below z=0.008 (although see § 3.4.3.). There is a 
marginally significant excess of absorber-galaxy pairs 
within 500 kpc, over that expected for a random dis- 
tribution of absorbers, although the inclusion of the 
Virgo velocity range removes the significance of this 
result. In summary, these tests seem to be consistent 
with the nearest-neighbour distance results, showing 
that there is an excess of close pairs of absorbers and 
galaxies, but that this result vanishes if the averaging 
is done over several galaxies or large radii. 

3.3.4- Tests of the 2nd Null Hypothesis: Identi- 
cal Lyman a absorber- Galaxy, Galaxy- Galaxy 
Correlations 

In order to test this hypothesis, one would like to 
use the same tests as were used in § p. 3. 3. A diffi- 
culty arises though in generating a large number of 
Monte Carlo samples. We are loathe to compare the 
observed distributions with simulations involving any- 
thing other than the actual observed galaxy distribu- 
tion since differences between the observations and 
simulations (based, for example on n-body or other 
galaxy clustering models) may result simply from in- 
adequacies of such models, and it is not clear how to 
create simulations of (fake) absorber- (real) galaxy dis- 
tributions having cross correlation properties which 
are the same as the observed galaxy-galaxy correla- 
tion properties. 

One way to deal with this difficulty is to use a test 
which does not require the generation of Monte Carlo 
samples: If the absorbers are distributed in the same 
way as the galaxies, then, given the pencil-beam na- 
ture of the galaxy sample, the redshift distributions 
of the absorbers and the galaxies should be identical, 
after correction for differing selection effects in the 
two samples. Table ^ shows the Kolmogorov-Smirnov 
D-values and probabilities that the absorber redshift 
distribution is the same as that of the galaxies, after 
the galaxy distribution is corrected using the selec- 
tion function shown in figure |3|. This selection func- 
tion was derived by assuming a Schecht cr luminos- 
ity funct ion with M*=-19.5 and a=-0.97 ( Loveday ei 
al. 1992 ). Because this test requires a known selec- 



tion function, it can only be run for the 'cone' galaxy 
sample. Also, as it directly compares the redshifts, it 
does not require any assumptions about pure or per- 
turbed Hubble flow. It can be seen from table ^ that 
when all the absorption lines are included there is a 
highly significant difference in the redshift distribu- 
tions. This significance level becomes marginal when 
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the Virgo absorbers are removed, and vanishes when 
only strong absorbers beyond Virgo are considered. 

Our other test of the hypothesis that the absorbers 
and galaxies have identical correlation functions uses 
the observed galaxy sample to generate our Monte 
Carlo 'absorber' sample. For each realization a num- 
ber of the actual galaxy redshifts were selected at 
random, and were treated as absorbers on the 3C273 
sightline. The actual algorithm involved selecting at 
random a number of the observed galaxy redshifts 
equal to the number in the absorber sample (making 
no correction for the galaxy selection function, and 
with no restriction on how close together the chosen 
galaxies were) , and treating these redshifts as if they 
were measured absorber redshifts on the 3C273 line 
of sight. The galaxies which provided these redshifts 
were removed from the galaxy sample for each test, 
to avoid an excess of 'spurious' pairs[| The number 
of absorber-galaxy pairs within a given radius were 
determined for both the real and 1000 Monte Carlo 
samples, in an identical manner to the second test in 
§ 3.3.3. above. This procedure should be valid as long 
as the radius within which the pair counts are being 
made is significantly larger than the typical distance 
between a galaxy and the 3C273 line-of-sight. Be- 
cause of this, tests were not run for a 500 kpc radius. 
The results from these tests are given in table 0. The 
columns are the same as for table |^, except that the 
final column lists the number of Monte Carlo runs 
with fewer absorber-galaxy pairs than the real sam- 
ple. Thus this number divided by 1000 is the prob- 
ability that the absorbers could have a correlation 
function as strong as that between galaxies. As can 
be seen, the absorber-galaxy correlation function av- 
eraged over 10 Mpc is significantly weaker than that 
between galaxies. 

In summary, the two tests above seem to show 
that the absorber-galaxy correlation function is sig- 
nificantly weaker than the galaxy-galaxy correlation 
function over large scales (10 Mpc). Even though 
there is clear evidence for galaxy-absorber clustering, 
there is a significant difference between the strength 
with which the galaxies are clustered with respect to 
each other and the strength with which the Lyman a 



absorbers are associated with galaxies. 

Some of the results of the two preceding sections 
can be inferred directly by inspection of the actual 
two-point correlation functions themselves. A log- 
arithmically binned version of the correlation func- 
tions for pure and perturbed Hubble flows are shown 
in figure [7| The absorber-galaxy correlation function 
was generated using the total Lyman a sample and 
the cone galaxy sample. The correlation functions 
were normalised in the usual way using random sam- 
ples with the selection function shown in figure [| 
The error bars were estimated using the formulae in 
Mo, Jing and Borner 1992. As expected, the pure 



e In fact leaving in the galaxies would allow one to model a situ- 
ation where the absorbers were actually part of the halo of one 
of the observed galaxies, although some maximum halo size 
would also have to be imposed to make the simulation realis- 
tic. We are investigating this point in morn detail and a paper 
is in preparation (Mo, Morris and Weymann 1993) 



and perturbed Hubble flow models agree fairly well 
for separations larger than one or two Mpc; the per- 
turbed Hubble flow model produces more very close 
pairs, since even small velocity differences wipe out 
small separations for the pure Hubble flow model. In 
both cases, the absorber-galaxy correlation is clearly 
weaker than the galaxy-galaxy correlation on scales 
from about 1-10 Mpc. However, although visual in- 
spection of the absorber-galaxy correlation function 
may suggest that the correlation is significant out to 
about 10 Mpc, in fact the pair tests summarized in 
table H and the nearest neighbor tests summarized in 
table H both indicate that a statistically significant 
absorber-galaxy correlations can be detected in our 
data sets only over volumes which are smaller than 
this. 

3.4. Some Particular Cases of Interest 

3.4.1. The Closest Absorber- Galaxy Associations 

Actual associations between individual observed 
galaxies and absorbers in the 3C273 sight line are dif- 
ficult to prove for a number of reasons: (1) The small- 
est projected distance to the 3C273 line-of-sight of all 
the galaxies in our sample is still 160 kpc. (2) Galaxy 
rotation or velocity dispersions could produce veloc- 
ity differences as large as 200 km/s between the mean 
galaxy velocity and an actually associated absorber 
(also comparable to the 3<r error in our fiber data 
velocity measurements) (3) In regions of high galaxy 
density along the line-of-sight, peculiar motions of the 
galaxies and the absorbers in cluster potential wells 
may make the velocity-distance relationship complex. 
This is especially true of the Virgo region, and for 
the z=0.0034 and 0.0053 absorbers (which also lie on 
the steep portion of the damping wings of the galactic 
Lyman a, and have only been observed at low resolu- 
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tion). With these caveats, it can be seen from tabic |] 
that the closest absorber-galaxy pair has a separa- 
tion of 350 kpc (240 kpc). Outside of the Virgo re- 
gion (where the velocity-distance relationship may be 
less complex, but also where our galaxy sample goes 
to much less faint absolute magnitudes), the smallest 
separation is 410 kpc (350 kpc). 

The best published example of an association be- 
tween a Lyman a absorber and a galaxy is given in 
Bahcall et al. 1992a, where they find a galaxy within 
90 kpc of the sightline to H1821+643 which has a 
strong Lyman a absorber (EW 950 mA) within 400 
km/s. There are no galaxies in our sample this close 
to the line of sight (or indeed any absorption systems 
this strong). 

One can (somewhat arbitrarily) break the absorbers 
in our sample into two groups: (1) those with a galaxy 
within 1 Mpc of the linc-of sight, and with a velocity 
difference of less than 400 km/s (all of the Virgo sys- 
tems and 6 of the highcr-rcdshift systems. Despite its 
entry in table |[ the absorber at z=0. 02622 actually 
has 2 galaxies within 400 km/s and 1 Mpc projected 
separation. They do not appear in the table, as their 
large velocity differences (>250 km/s) make their sep- 
aration large, assuming pure Hubble flow), and (2) the 
rest (8 systems). There is no significant difference in 
the EW distribution of these two samples. 

Our data can also be used to consider the ques- 
tion: What is the average galaxy diameter within 
which one would see a neutral-hydrogen column den- 
sity of at least 10 13 cm~ 2 ? In practice this is rather 
a naive question, as the cross-section almost certainly 
depends on galaxy luminosity and probably also mor- 
phology. One might also expect a patchy distribu- 
tion of neutral-hydrogen in the outer parts of galaxies 
leading to a covering factor not equal to one. Never- 
theless, ignoring these complications, and excluding 
the Virgo velocity range due to (a) the possibility of 
large peculiar velocities (b) a wish to avoid the large 
range of intrinsic galaxy luminosities in the sample 
and (c) because our absorber EW limit is higher in 
this region, we find that of the 12 galaxies with pro- 
jected separation to the line-of-sight less than 1 Mpc, 
8 show Lyman absorption systems within 400 km/s. 
A velocity difference of up to about 400 km/s could 
possibly be attributed to internal motions within a 
large galaxy, coupled with our measuring error, or al- 
ternatively to a small group. Having 8 or more such 
matches would occur 0.6% of the time if the absorbers 
were randomly distributed in velocity space. How- 



ever, of these 8 galaxies, two are associated with the 
absorber at z=0. 02622, and three with the absorber 
at z=0.2933. Clearly several distinct galaxies can not 
be producing the same absorber, and so in fact only 5 
of the 12 systems within 1 Mpc of the line-of-sight to 
the quasar can be legitimately associated with indi- 
vidual absorbers. Reducing the projected separation 
to 500 kpc, one finds 3 galaxies, with velocity sep- 
arations of 20, 120 and km/s (see table ||). This 
may be interpreted as saying that the cross-section of 
galaxies with luminosities greater than 1/10 of L* for 
log(NH)> 10 13 is between 0.5 and 1 Mpc, although 
by "cross-section" we do not mean that the absorbers 
are neccesarily associated with the actual galaxy in 
question. 

Because of the 3 Virgo absorbers, all of the Virgo 
galaxies have an absorber within 400 km/s. However, 
it is clearly unreasonable to suggest that more than 
one galaxy is associated with a given absorber. Thus 
of the 19 Virgo galaxies within 500 kpc of the line- 
of-sight, 16 must not be producing absorption. Of 
these 19, 2 have absolute magnitudes brighter than 
1/10 L*, and so are consistent with the above state- 
ment, despite the higher absorber EW limit, as long 
as only one of the galaxies with luminosity < 1/ 10 L* 
is producing observed absorption. 

3.^.2. The Most Isolated Lyman a Absorbers 

One of the stronger lines in our absorber sample 
is also the most isolated. Any galaxy brighter than 
-18 would be in our sample at the redshift of the iso- 
lated z=0. 06655 absorber. However, the nearest such 
galaxy has a projected separation of 4.75 Mpc and 
a velocity difference of 710 km/s, which corresponds 
to a spatial separation of 10 Mpc for both pure and 
perturbed Hubble flow models. Indeed studying the 
galaxy distribution in figure ^ the three absorbers 
with 0.060<z<0.07 seem to lie in a 'void' in the galaxy 
distribution. While it is certainly possible that these 
absorbers are associated with galaxies below our ab- 
solute magnitude limit, these isolated absorbers are a 
clear demonstration that one can not associate ALL 
low-redshift Lyman a absorbers with luminous (L* ) 
galaxies. 

3.4-3. The Absence of Lyman a Absorbers in the z 
~ 0.078 Galaxy Concentration 

The most striking feature of the observed galaxy 
redshift distribution is a concentration of galaxies cen- 



13 



tered at a redshift of z ~ 0.078. This excess is almost 
certainly associated with a structure which includes 
the galaxy cluster Abell 1564. This cluster has a tab- 
ulated center at 12:32:25 +02:07:11 (1950), giving it 
an offset of 89' SW of the 3C273 sightline. This angle 
corresponds to a separation of 7.3 Mpc at the clus- 
ter distance. Abell 1564 is the closest Abell cluster 
to the 3C273 sightline, but is only richness class 0. 



Redshifts for two member galaxies are given in Met 



cale et al. 1989, giving an average cluster redshift of 



0.0793. Selecting all galaxies in our fiber sample with 
redshifts between 0.07 and 0.085 (54 galaxies), one 
derives a mean redshift of 0.0781 with rms of 0.0021 
(630 km/s). For comparison, selecting galaxies be- 
tween redshifts of 0.15 and 0.17 (the cluster around 
3C273, see figure ||), one gets 24 galaxies with mean 
redshift 0.1581 and rms 0.0024 (720 km/s). 

Despite the location of Abell 1564 to the SW, the 
galaxies in the redshift slice 0.07-0.085 show a weak 
concentration to the NE of the fiber field, although 
there are galaxies in all parts of the fiber area sur- 
veyed. If they are really part of a structure includ- 
ing A1564, the structure must have a size of at least 
8 Mpc. However, the nearest Lyman a absorber is 
0.0095 from the peak in redshift space (or 2850 km/s, 
which is 4.5 times the rms from the mean of the con- 
centration) and is at a distance of 36 Mpc from the 
center of the concentration. It thus seems unlikely 
that there are any observed Lyman a absorbers in 
our line of sight which are physically associated with 
this concentration. 

We can test a slightly different form of the hypoth- 
esis considered in § |3.3.4. , (though we grant its appli- 
cation to the z ~ 0.078 galaxy concentration involves 
post facto statistics), namely: What is the probabil- 
ity that, if the Lyman a absorbers are distributed in 
space "in the same way" as the galaxies in our sam- 
ple, we should find the observed number of Lyman a 
absorbers (none, in our case). By "in the same way" 
we mean that: 



PLy a (r) 



Phyu 



P galaxy 

if) 
Pgalaxy 



and by "in the vicinity of" we shall mean within 
±2.5 x rms of the redshift peak of the concentration. 
We consider that our selection function can be mean- 
ingfully applied over the redshift range from about 
z= 0.016 out to our adopted "proximity cutoff" at 
z= 0.151. Carrying out the integration of our sam- 
ple (i.e. the galaxy numbers weighted by the selection 



function) between ±2.5 x rms of the peak of the galaxy 
concentration and over the redshift range from 0.016 
to 0.151 we find that after correction for the selection 
function, about 0.33 of an absolute-magnitude-limited 
cylindrical sample of galaxies should be found within 
2.5cr of the velocity peak of this structure. Since the 
Lyman a absorbers are not subject to such a selec- 
tion function, we can apply this equation to all the 
volume elements along the line of sight and, assuming 
further that the effective cross-section of the Lyman 
a absorbers (<r) is constant over the redshift range 
0.016-0.151, one can write: 

n x a x [(z peak + 2.5 x rms) - (z. peak - 2.5 x rms)] 

= 0.33 x WTt x [0.151 - 0.016] 

This equation represents the number of Lyman a ab- 
sorbers in the vicinity of the peak that we should ex- 
pect to see if our hypothesis is correct. Over the red- 
shift interval 0.016-0.151 we have observed up to 14 
Lyman a absorbers so that we expect 0.33 x 14 as 4.7 
absorbers, and we observe none. The probability of 
this occuring is thus exp(-A.T) w 0.01. This is not a 
strong result due to the rather small number of ab- 
sorbers expected, and the post facto nature of the 
test as we acknowledged above, but is additional sug- 
gestive evidence that the Lyman a absorbers do not 
follow the galaxy distribution and in particular that 
they may avoid strong concentrations of galaxies. It 
could be strengthened, of course, by observing addi- 
tional lines of sight through dense concentrations. 

It would be nice to repeat the above calculation 
for the case of the three absorbers found in the Virgo 
velocity region. However, there are a number of rea- 
sons why any estimate of the probability of finding 
three absorbers in this region with an apparent over- 
density of galaxies is highly uncertain. Firstly, our 
absorber sample is inco mplete in this velocity range, 
as described in 



3.3.1 



Secondly, the area of sky sur- 
veyed is so small that the selection function correc- 
tion for this velocity range is very large. Taking our 
selection function at face value, the 8 galaxies actu- 
ally in our survey at Virgo velocities imply that 60% 
of an absolute-magnitude-limited cylindrical galaxy 
sample from z=0 to 0.151 would be found in that 
velocity range. Our hypothesis above would then 
predict that (neglecting the incompleteness of the 
absorber sample) one would expect 0.6xl7«10 ab- 
sorbers, when we only see 3. Unfortunately, it is not 
possible to place a meaningful uncertainty on this ap- 
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parent under-density of absorbers in the Virgo region, 
for the reasons listed above. 

In summary, there is marginal evidence, both from 
a clump of galaxies at z=0.078 and also from the 
Virgo velocity range, that Lyman a absorbers are less 
common in regions of high galaxy density. 

4. Summary and Discussion 

We have assembled several types of observations in 
an attempt to find objects with which the Lyman a 
absorbers along the line of sight to 3C273 might be 
associated, and in order to carry out statistical tests 
of galaxy-Lyman a absorber association. In particu- 
lar, we obtained narrow-band images centered on and 
off the expected position of any Ha emission which 
might be associated with three of the low redshift Ly- 
man a absorbers, and obtained deep broad-band im- 
ages of a 17'xl7' field centered on 3C273. Both these 
searches were negative. Our failure to identify any 
broad-band or Ha emission from plausible "galaxy- 
like" objects a few tens of kpc from the 3C273 sight 
line at the approximate distance of the Virgo clus- 
ter will be checked by more sensitive and extensive 
searches for Ha by T. Williams (private communica- 
tion) and 21 cm emission by |van Gorkom 1993 and 
van Gorkom et al. 1993. 



We have also obtained redshifts for a large num- 
ber of galaxies in the vicinity of the 3C273 sightline. 
Again, we find no unambiguous instance of associa- 
tion of any of the Lyman a absorbers with individual 
galaxies. We define a number of samples for both the 
Lyman a absorbers and the galaxies and estimate the 
3-dimensional separation between each galaxy-galaxy 
pair and each Lyman a absorber-galaxy pair based 
upon two models for converting the observed redshift 
difference between any pair into a radial separation, 
viz. i) the assumption of a pure Hubble flow and ii) 
a statistical m odel of "perturbed" Hubb le flow based 
upon work of Davis and Peebles 1983. The result- 



ing data base is used to carry out statistical tests to 
confirm or reject two null hypotheses about the asso- 
ciation of galaxies and Lyman a absorbers, namely: 
i) The Lyman a absorbers show no tendency to clus- 
ter around galaxies ii) The Lyman a absorbers cluster 
around galaxies exactly as the galaxies cluster about 
each other. While neither of these two hypotheses 
can be unambiguously rejected in the sense that ev- 
ery combination of samples and flow hypotheses reject 
both of them at significant levels, the evidence from 



these tests, and from the galaxy-galaxy and Lyman 
a absorber-galaxy two-point correlations themselves, 
points quite strongly to the conclusion that both hy- 
potheses are false. 

In particular, over length scales from about 1 to 10 
Mpc there seems little doubt that the Lyman a ab- 
sorbers cluster around galaxies less strongly than the 
galaxies themselves cluster. This is born out by an 
examination of a redshift interval centered at about 
z=0.078 at which a strong concentration of galaxies 
occurs but in the neighborhood of which there are 
no Lyman a absorbers. Additionally, we find at least 
one Lyman a absorber for which no galaxy with abso- 
lute magnitude brighter than about -18 can be found 
closer than about 5 Mpc. Taken together, all this ev- 
idence suggests that the most significant conclusion 
we have reached is that the majority of low-redshift 
Lyman a absorbers are not intimately associated with 
normal luminous galaxies. 

In view of the fact that it has long been realized 
that at high redshifts there is an absence of power in 
the Lyman a absorber two-point correlation function 
in redshift space, except possibly at the very smallest 
velocity separations, this conclusion is not too sur- 
prising. On the other hand, it is also fairly clear that 
there is some tendency for the Lyman a absorbers to 
cluster around galaxies, and even weak evidence that 
this clustering becomes strong at very small separa- 
tions. Also Bahcall et al. 1992a have investigated 



the auto-correlation function of the low-redshift ab- 
sorbers seen in the line-of-sight to H1821+643, and 
show that there is only a 4% probability that the ob- 
served 'clumping' arose from a randomly distributed 
sample. 

None of the foregoing points unambiguously, in our 
estimation, to a particular model for the formation 
and evolution of the Lyman a absorbers. As have oth- 
ers, we simply offer the following speculations which 
appear to be compatible with the facts as they are 
presently understood: 

At high redshifts, the Lyman a absorbers consist 
primarily of entities which are only very loosely asso- 
ciated with larger mass objects (e.g. proto galaxies), 
and which are evolving fairly rapidly. Possibly this 
dominant population consists of absorbers in which 
gravitational binding (eg by dark matter) plays no 
significant role. In addition to this group, there is 
a smaller population of absorbers which are evolving 
less rapidly, possibly stabilized by dark matter and 
which are clustered more strongly about galaxies. At 
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very low redshifts, this latter population is beginning 
to constitute a large enough fraction of the absorbers 
that power in the two-point correlation function is 
detectable. Thus, the present mix of Lyman a ab- 
sorbers appears to have clustering properties inter- 
mediate between present-epoch normal galaxies and a 
random non-clustered population. This property also 
appears to be shared by the low-luminosity moderate 



redshift "blue galaxies" (Pritchct and Infante 1992), 
leading to the plausible conjecture that the Lyman a 
absorbers are more closely related to low mass, low lu- 
minosity galaxies than they are to L* galaxies, though 
the relation is clearly not one-to-one. 

We have no good way at present of estimating the 
characteristic scale or masses of the low-redshift Ly- 
man a absorbers. A guess at a diameter of 30 kpc 
is as plausible as any. In particular, consider a pan- 
cake whose diameter is 30 kpc and whose thickness is 
10 kpc. In this case, at the present epoch, a hydro- 
gen column-density of 10 13 -10 14 cm~ 2 normal to the 
face of the absorber, coupled with estimates for the 
present-epoch energy-density of ionizating radiation 
leads to a total gas mass of order 10 7 M , but a mass 
of only a few hundred solar masses of neutral hydro- 
gen. It is of interest that the mass function of neutral- 
hydrogen gas clouds appears to be truncated below 
about 10 s Mq flWeinberg et al. 199l| ). As Maloney 
(1992,1993) has shown, the gas in a flaring galaxy 
with decreasing column-density will undergo a rather 
sudden transition along its face from being mostly 
neutral to mostly ionized, with the consequence that 
few if any contours with neutral-hydrogen column- 
density of order 10 18 cm -2 are known. Similarly, it is 
conceivable that, given the appropriate run of length 
scale with mass, a sequence of masses would have 
the property of making a sudden transition in the 
mass function of neutral-hydrogen starting at about 
10 8 Mq, leading to a dearth of objects with total HI 
masses for several orders of magnitude below this. 

If these speculations have any connection with re- 
ality then one might expect to see some similarity in 
the clustering properties of the Lyman a absorbers 
and low mass galaxies. We are currently attempt- 
ing to obtain redshifts of galaxies of lower luminosity 
along the 3C273 sightline in order to investigate this 
possibility. 
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Table 1: Definite Galaxy redshifts in 3C273 Field 
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Table 2: Comparison of 3C273 Extragalactic Absorption Line Lists 
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Table 3: Nearest Neighbour Monte-Carlo tests 



Table 4: Nearest Galaxies to 3C273 Lyman a absorbers 
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Table 5: Pair Count Monte-Carlo tests using random 



Table 6: KS comparison of z- distribution 



Table 7: Pair Count Monte-Carlo tests using the galaxy distribution 
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Fig. 1. — Coronagraph Data of 3C273 field, (a) Sum 
of 25A Filters, the horizontal bar in the lower left is 1 
arminlong, (b) Continuum Subtracted VN1, (c) Con- 
tinuum Subtracte d VN 2, (d) Continuum Subtracted 



HN. See text in 5 2.1 



Fig. 2. — Mosaic of CO SMI C Images in Gunn r of 
3C273 field. See text in § 2/2. Region shown in figure [l] 



is circled. 



Fig. 3. — Histogram of redshifts for Galaxies found in 
the 3C273 sightline. Also plotted are the locations of 
Lyman a absorbers from table || (where 'strong' refers 
to lines with EW>55 mA, and 'weak' refers to the 
remainder), and the (arbitrarily normalised) galaxy 
selection function for the fiber survey. See text in 



2.3 



Fig. 4. — Locations on the sky of the galaxies with 



redshifts in the field of 3C273. See text in 8 2.3 



Fig. 5. — Pie-diagrams for the Galaxies observed with 
the LCO fiber system. Angles have been exagerated 
by a factor 15 to prevent overcrowding of the symbols. 
Please note that this results in a highly distorted plot 
with initially spherical structures (such as the 3C273 
cluster of galaxies) appearing elongated transverse to 
the line of sight. Note also that the star marking the 
position of 3C273, while readily visible in the pro- 
jection in RA, is partially obscured by a clump of 
galaxies in the projection in Dec 



Fig. 6. — Mosaic of COSMIC Images in Gunn r of 
3C273 field, with 3 artificial dwarf galaxies added. 
NE with M B =-14.5 and scale length 15kpc, NW with 
Mb=-13.5 and scale length 7.5kpc, and SW with 



Mb=-13.5 and scale length 15kpc. See text in § 3.2. 



Fig. 7. — Two point correlation functions for galaxy- 
galaxy and absorber-galaxy (a) assuming pure Hubble 
flow (b) assuming perturbed Hubble flow. See text in 



3.3.4 
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